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ABSTRACT 


Several  different  statistical  fracture  theories  are  developed  for 
materials  with  cracks  confined  to  the  surface.  All  assume  that  crack  planes 
are  normal  to  the  surface,  but  are  otherwise  randomly  oriented.  The  simplest 
theory  assumes  that  only  the  component  of  stress  normal  to  the  crack  plane 
contributes  to  fracture.  This  theory  is  in  fair  agreement  with  biaxial  frac- 
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ture  data  on  pyrex  glass  obtained  by  Oh.  When  the  contribution  of  shear  is 
included  in  the  analysis,  the  crack  shape  has  to  be  considered.  Several  shapes 


are  examined,  and  the  corresponding  fracture  statistics  are  derived 


ure  criterion  employed  is  that  fracture  occurs  when  the  maximum  tensile  stress 


rial.  The  assumption  of  shear-sensitive  cracks  leads  to  improved  agreement 
with  experiment,  but  really  good  agreement  appears  to  require  the  assumption 
that  the  cracks  have  a preferred  orientation. 
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Fig.  1 Probability  of  failure  for  pyrex  tubes  under  biaxial  tensile 
stress  states  1:1,  1:0.5,  and  1:0  assuming  only  normal  com- 
ponents of  stress  on  the  cracks  contribute  to  failure.  The 
curves  for  1:0.5,  and  1:0  are  generated  from  the  curve  for 
1:1,  which  is  fitted  to  the  experimental  results  of  Oh  (1970). 

Data  points  for  the  three  stress  states  are  plotted 6 


Fig.  2 Failure  probability  curves  fitted  to  1:1  data  of  Oh  (1970). 

Curve  (0)  is  determined  from  parameters  reported  by  Oh  for  the 
three-parameter  Welbull  function  Pf  ■ 1 - exp{-  ((a  - CTU) /ctq)™} . 
Curve  (A)  is  a Ueibull  function,  and  (B)  is  a piecewise  function 
consisting  of  two  Welbull  curves  Joined  by  a straight  line  such 
that  Pf  and  its  slope  remain  continuous.  The  values  of  the 
parameters  are  as  follows: 
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Fig.  3 Probability  of  failure  for  uniaxial  tension  by  Oh's  method. 

The  curves  labeled  A(1:0),  B(1:0)  and  0(1:0)  are  Welbull 
functions  whose  parameters  are  determined  graphically  from 
computed  values  of  Pf  « 0.01  starting  from  the  A(l:l), 

B(l:l),  and  0(1:1)  fits  to  the  equlblaxlal  data,  respec- 
tively. The  unlabeled  curve  was  plotted  using  the 
Welbull  parameters  for  1:0  reported  by  Oh. 10 


Surface  crack  models:  (a)  Griffith  through-crack, 

(b)  Griffith  notch,  (c)  half  ellipsoid.  Pj_,  ?2,  P3  indi- 
cate the  locations  of  the  maximum  stress  in  each  model. 

P3  lies  in  the  X],  - X3  plane,  but  not  on  the  X^  axis  when 
shear  is  present.  Its  exact  position  depends  on  the 
stress  state  


Failure  probability  curves  for  uniaxial  tension  (1:0) 
determined  from  the  B(l:l)  fit  to  the  equlblaxlal  data 
assuming  three  different  crack  models:  the  shear 

insensitive  (SI),  the  Griffith  notch  (GN) , and  the 
half-ellipsoid  (HE).  The  points  are  the  1:0  data  of 
Oh  (1970) 
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I.  INTRODUCTION 


There  are  many  cases  in  which  structures  must  be  fabricated  using  brittle 
materials.  The  strongest  and  most  refractory  materials  tend  to  be  brittle. 
Also,  materials  transparent  to  microwave,  infrared,  or  visible  radiation  are 
generally  brittle.  Brittle  structures  characteristically  exhibit  a large 
variation  in  fracture  stress  which  must  be  taken  into  account  in  design. 

The  most  widely  used  statistical  theory  of  fracture  is  due  to  Weibull 
(1939).  He  attributed  the  variation  in  fracture  stress  of  nominally  identical 
specimens  to  the  presence  of  unidentified,  invisible  flaws.  The  flaws  were 
assumed  to  have  a distribution  in  strength,  and  the  specimen  or  structure  was 
assumed  to  fail  when  the  strength  of  the  weakest  flaw  or  link  was  exceeded. 

Batdorf  and  Crose  (1974)  revised  weakest  link  theory  by  assuming  the 
flaws  to  be  cracks,  and  therefore  to  have  strengths  which  depend  on  the  orien- 
tation of  the  cracks  with  respect  to  the  applied  stresses.  All  orientations 
were  considered  equally  likely,  i.e.,  the  material  was  assumed  to  be  macro- 
scoplcally  Isotropic.  It  was  further  assumed  that  only  the  component  of 
stress  normal  to  the  crack  plane  contributed  to  fracture.  The  latter  assump- 
tion was  a convenient  approximation  which  permitted  development  of  a general 
theory  without  having  to  specify  crack  shapes.  The  shear  parallel  to  the 
crack  plane  also  contributes  to  the  fracture,  but  by  an  amount  that  depends 
on  the  crack  shape,  which  is  something  one  usually  does  not  know.  In  some 
cases,  however,  we  may  be  able  to  derive  information  about  crack  shapes  by 
examining  the  fracture  statistics  for  a number  of  different  stress  states. 

Only  volume  distributed  cracks  were  treated  by  Batdorf  and  Crose  (1974). 
In  the  case  of  some  materials,  e.g.,  glass.  It  is  generally  accepted  that 
*11  cracks  are  located  at  the  surface,  and  also  that  the  crack  planes  are 
normal  to  the  glass  surface  (McCllntock  and  Argon,  1966) . For  such  materials 
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the  crack  orientation  is  given  by  a single  parameter  rather  than  two,  as  in 
the  case  of  volume  distributed  cracks. 

In  the  present  paper  the  theory  of  Batdorf  and  Crose  is  modified  for 
surface  distributed  cracks  and  is  applied  to  fracture  data  for  glass  obtained 
by  Oh  (Oh,  1970;  Oh  et  al.,  1973).  The  agreement  with  experiment  leaves  some- 
thing to  be  desired.  It  is  evident  that  Including  the  effects  of  shear  on 
the  crack  plane  would  decrease  the  discrepancy. 

A more  refined  theory  including  the  Influence  of  shear  and  based  on  the 
assumption  that  the  cracks  are  Griffith  cracks  and  that  failure  occurs  when 
the  local  tensile  stress  on  the  crack  surface  exceeds  the  intrinsic  strength 
of  the  material  was  developed  by  Oh  (1970).  Some  improvement  in  agreement 
with  test  data  resulted.  The  present  authors  believe  that  the  crack  model 
employed  by  Oh  is  not  appropriate  for  surface  cracks.  Alternative  models 
are  therefore  proposed,  but  the  improvement  is  marginal.  Good  agreement 
with  experiment  is  obtained  by  assuming  that  in  addition  to  being  shear- 
sensitive,  the  cracks  have  a preferred  orientation. 


II.  THEORY:  SHEAR  ON  CRACK  PLANE  NEGLECTED 

* 

Consider  first  a single  crack  of  arbitrary  orientation.  Since  it  is 
assumed  to  be  small  and  located  at  the  surface,  it  is  subjected  at  most  to 
plane  stress.  In  the  principal  axis  system,  the  tensile  component  of  stress 
normal  to  the  crack  line  and  in  the  plane  of  the  surface  is 
o 


2 2 
a cos  0 + a sin  0 
n x y 


(1) 


where  0 is  the  angle  between  the  x-axis  and  the  crack  normal.  In  accordance 


with  the  preceding  assumptions,  the  material  will  rupture  when  or  > ocr 

where  ocr  is  the  macroscopic  normal  stress  required  to  rupture  the  crack. 

If  the  crack  is  randomly  oriented,  the  probability  of  failure  is  given  by 

P,  - « 

f ir 


(2) 
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where  w is  the  radian  measure  of  the  angular  range  in  the  positive  ox  half- 
plane within  which  o > a . 

n cr 

In  a real  material,  there  will  be  a number  of  cracks  of  varying  orienta- 
tion and  critical  stress.  If  we  assume  that  the  cracks  are  uniformly  dis- 
tributed over  the  surface,  the  material  can  be  characterized  by  a density 
function  °cr  ) where  is  the  applied  stress  state.  This  function 

represents  the  number  of  cracks  per  unit  area  having  a critical  stress  less 
than  or  equal  to  Ocr<  The  number  of  cracks  per  unit  area  having  critical 

stresses  between  o and  o + da  is,  then, 
cr  cr  cr 


dN 


dN 

-j- — do 
dacr  cr 


(3) 


The  probability  that  failure  will  occur  in  a uniformly  stressed  surface  of 

area  A due  to  a crack  having  a critical  stress  in  the  range  a to  a + da 

cr  cr  cr 

is  the  product  of  the  probability  that  a crack  is  present  and  the  probability 
that  the  crack,  if  present,  will  fail;  i.e., 

Pf(E  ■ d0cr)  - (‘  df;  ■“’cr)  (if) 

The  probability  that  such  cracks  will  survive  is 

P.(£.  ‘“’cr)  “ 1 “ p 

a exp 


(A) 


_ 1 «Ld0 

f tt  do  cr 

cr 


r_A£^N_da  i 

L 77  d°cr  crJ 


(5) 


The  probability  that  cracks  in  every  stress  range  da  will  survive  is  the 

cr 

product  of  the  probabilities  of  survival  of  cracks  in  the  individual  ranges, 
i.e. , 

r t M HN  1 

p. 


Pg  - exp 


[-*/?#;  docr]-  »' 


(6) 
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The  fracture  probability  for  a surface  of  area  A subjected  to  stress 

state  ^ can  be  evaluated  with  the  use  of  Eq.  (6)  when  to  and  N are  known. 

For  any  stress  state  , we  can  determine  U)  by  using  Eq.  (1),  while  N must 

be  obtained  by  experiment. 

In  the  uniaxial  case,  Eq.  (1)  reduces  to 

2 

o - o cos  9 
n x 

Thus  for  this  case,  ±9cr»  the  angle  within  which  On  > Ocr>  is  given  by 

2 

o - a cos  9 
cr  x cr 


or 


Since  u)  * 20 


0 * cos 

cr 


cr’ 


u>  2 -1  / gcr 

— * — COS  \ / 

7T  7T  V O 


In  the  equlbiaxial  case. 


0*0-0 
x y 


as  a result  of  which,  o - o and 

n 


u> 


£ - 1 for  o <o 
it  cr 


- 0 for  o >0 
cr 


(7) 


(8) 

(9) 

(10a) 

(10b) 


As  a result  of  (10a,  b),  in  the  equlbiaxial  tension  case  Eq.  (6)  takes  the 
simple  form 

Pf(o)  - 1 - exp  [ - AN(o) ] (11) 
From  a comparison  of  Eqs.  (6)  and  (8)  with  (11),  it  becomes  clear  that  it  is 
computationally  advantageous  to  determine  N from  the  fracture  statistics  for 
equal  biaxial  tension.  From  (11),  we  find  AN  is  given  by 


J 


£ 


(12) 


thereby  avoiding  the  necessity  encountered  in  the  theory  of  volume-distributed 
cracks  of  solving  simultaneous  linear  algebraic  equations  or  an  Integral  equa- 
tion to  obtain  N.  Actually,  in  the  volume  distribution  case,  a simplifica- 
tion like  Eq.  (12)  occurs  for  equltriaxlal  tension,  but  this  state  of  stress 
cannot  be  realized  in  practice. 

Under  general  biaxial  stress  states,  Eq.  (1)  can  be  rewritten  as 

°n  " ax(CO820  + K sin29)  * ax[cos20(l  - K)  + K]  (13) 

where 


o 


(15) 

(16) 


Using  Eqs.  (6),  (12),  and  (16)  we  can  obtain  the  probability  of  failure  of  a 


surface  of  area  A subject  to  biaxial  tension. 


The  results  of  the  theory  just  outlined  are  compared  with  the  experi- 
mental results  of  Oh  in  Figure  1.  It  is  evident  that  although  the  general 
trend  of  the  theoretical  results  are  in  accord  with  the  data,  the  theoretical 
curves  for  stress  ratios  1:1,  1:0.5,  and  1:0  are  too  far  apart.  Taking 
account  of  the  contribution  of  shear  on  a crack  to  the  failure  process  would 
increase  the  probability  of  failure  for  stress  ratios  1:0.5  and  1:0,  and 
therefore  bring  the  curves  closer  together.  He  therefore  turn  to  theories 
in  which  shear  effects  are  taken  into  account. 
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III.  OH'S  THEORY 

It  was  shown  by  Oh  (1970)  that  when  a Griffith  crack  is  subjected  to 

principal  stresses  0 and  0 , the  maximum  tensile  stress  at  any  point  on  the 

x y 

surface  of  a crack  is  (in  our  notation) 


- £ +/r2  + S2  ^ 


where  the  tensile  and  shear  forces  on  the  crack  are  given  by 


(17a) 


(17b) 


S - 0 


sin  26 


(17c) 


In  these  equations  K = 0 /o  < 1,  £ is  the  ratio  of  the  minor  to  the  major 

y * 

axis  of  the  ellipse,  while  6 is  the  complement  of  the  angle  between  the  larger 
principal  stress  and  the  crack  plane.  All  values  of  6 were  assumed  equally 
likely,  and  the  distribution  of  cracks  with  respect  to  £ was  chosen  to  fit  a 
three-parameter  Weibull  representation  of  the  data  for  the  stress  ratio  1:1. 

To  obtain  the  failure  probabilities  for  stress  ratios  1:0.5  and  1:0  it 
was  arbitrarily  assumed  that  these  would  also  be  three-parameter  Weibull  dis- 
tributions. The  corresponding  parameters  were  found  by  applying  the  graphical 
approach  for  Weibull  parameter  estimation  to  synthetic  data  computed  for 
Pf  « 0.01.  Such  a procedure  would  be  theoretically  justified  if  failure 
always  occurred  for  a condition  not  applying  to  Oh's  data.  In 

spite  of  this,  when  the  procedure  is  applied  using  Oh's  Weibull  parameters  for 
equlbiaxlal  tension,  good  agreement  is  obtained  with  uniaxial  test  data.  How- 


ever, two  other  fits  to  the  equlbiaxlal  data  were  devised  by  the  present 
authors  (Figure  2).  One,  labeled  A,  la  an  alternative  Weibull  function.  The 
other,  labeled  B and  Identified  In  Figure  2,  Is  a closer  fit  to  the  data,  but 
also  more  complicated  than  any  Weibull  function.  Although  the  three  fits  to 
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Failure  Probability  Curroa  Finad  to  1:1  Data  of  Oti  (1970).  Curw  (0)  it  datarmtoiad  from 
paramatart  reportad  by  Oft  for  tha  thraa-paramatar  Waibull  function  Pf  “ 1 - axpj - ((a-ou)/oo)m } 
Curre  (A)  ia  a Waibutl  function,  and  (B)  to  a piacawiu  function  comtoting  of  two  Waibull  curvai 
toinad  by  a »twi#it  lina  tueh  that  Pf  and  its  dopa  remain  continuous.  The  valuat  of  the  poramotart 


the  1:1  data  are  quite  close,  the  Oh  procedure  leads  to  very  significant  dif- 
ferences In  the  predicted  uniaxial  failure  curves  for  the  three  cases,  as 
shown  In  Figure  3.  This  suggests  that  the  good  agreement  found  using  Oh's 
function  may  be  fortuitous. 

Oh's  analysis  Is  based  on  the  tacit  assumption  that  the  Griffith  crack 
is  a through-crack  (Figure  4a) . This  is  not  considered  by  the  present  authors 
to  be  an  appropriate  model  for  the  surface  crack.  Section  IV  is  concerned 
with  two  other  types  of  cracks  illustrated  in  Figures  4b  and  4c.  The  first, 
which  we  shall  call  a Griffith  notch,  is  a half-elliptic  cylinder  with  the 
principal  axis  of  the  cylinder  at  the  surface  of  the  specimen.  The  second  is 
a half -ellipsoid  in  which  a » b » c. 

IV.  PRESENT  THEORY 


A fracture  criterion  under  combined  stress  can  be  formulated  on  any  of 
several  different  bases.  These  include  use  of  an  energy  criterion,  use  of 
critical  stress  concentration  factors,  and  use  of  maximum  tensile  stress 
occurring  at  a point  on  a surface  of  the  crack.  The  simplest  one  of  these 
to  apply  is  the  last. 

Mlrandy  and  Paul  (1975;  also  Paul  and  Mirandy,  1975)  have  recently 
worked  out  the  stress  state  at  any  point  on  the  surface  of  an  ellipsoidal 
cavity  having  axes  a,  b,  and  c,  such  that  c « b £ a for  arbitrary  applied 
stresses.  For  present  purposes  the  principal  findings  are:  (1)  If  the 

applied  stress  is  simple  tension  T normal  to  the  plane  of  the  crack  (l.e., 
parallel  to  the  polar  or  c-axls  of  the  crack) , the  maximum  stress  occurs  at 
the  intersection  of  the  cavity  and  the  a-b  or  equatorial  plane.  The  local 
stress  is  the  same  at  all  points  in  the  equatorial  plane  and  is  given  by 


0 


b 2T 
c E 


(18) 
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2 
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«•«»  3.  Probability  of  Fadura  for  Uniaxial  Tanaion  by  Oh'i  Mathod.  Tha  curvas  labatad  A(1:0) 
B(1:0)  and  0(1:0)  ara  WaMI  functions  whosa  pammators  art  datarminad  graphically 
from  eomputad  mluaa  of  Ff4C0.01  ttartint  from  tha  A(1:1),  8(1:1),  and  0(1:1)  fits  to 
tha  aquibiaxial  data,  raspacthrsfy.  Tha  unlabaiad  ourva  sms  ptottad  using  tba  Waibutl 
paramatars  f or  1:0  raportad  by  Oh. 


Figure  4.  Surface  Crack  Models:  (a)  Griffith 
Through -Crack,  (b)  Griffith  Notch, 
(e)  Half  ENipaoM.  Pi,  P2.  P3  indicate 
the  locatiom  of  the  mexirnuui  atteat 
in  each  model.  P3  liea  in  the  X1-X3 
ptene,  but  not  on  the  Xi  axia  when 
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where  E is  an  elliptic  integral  of  the  second  kind.  For  b/a  « 1,  E - 1. 

(2)  If  shear  stress  S is  applied  to  the  crack  plane  in  a direction  parallel 
to  the  a-axis,  the  maximum  tensile  stress  induced  on  the  surface  of  the 
cavity  occurs  in  the  a-c  plane,  and  is  given  by 

(19) 


_ _ b 2 sin  26  _ 

° c E 2(1  - V)  S 


where  $ is  the  local  latitude.  (3)  When  both  S and  T as  described  above  are 
present,  the  local  tensile  stress  on  the  cavity  surface  is  largest  in  the  a-c 
plane  and  becomes,  for  a » b. 
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T cos  26  - 


2(1  - v) 


sin  2g 


which  has  a maximum  value  of 


max 


(1  - v)U 


(20) 


(21) 


(4)  In  the  case  of  a Griffith  crack  subjected  to  tension  normal  to  the  crack 
plane  and  shear  on  the  crack  plane  applied  parallel  to  the  cylinder  axis,  the 
tensile  stress  under  combined  loads  is  greatest  along  the  line  of  maximum 
curvature.  It  is  evaluated  by  finding  the  maximum  principal  stress  at  the 
end  of  the  b-axis  in  a long  ellipse,  where,  according  to  Hirandy  and  Paul 
(1975),  the  local  stress  state  la 
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(22a) 

z c 

o • v o (plane  strain) 

x z 

(22b) 

T - - S 

(22c) 

c 

Using  Mohr's  circle,  this  is  readily  shown  to  be 

r *l 

°max  " c lT(1  + V)  ‘•Vt2*1  * v>2  + S2  J 

(23) 
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The  above  results  are  for  ellipsoids  in  an  isotropic  elastic  body 
located  far  from  any  free  surfaces.  In  the  case  of  surface  cracks  such  as 
those  shown  in  Figure  4,  the  presence  of  the  free  surface  will  modify  to 
some  degree  the  stresses  on  the  surface  of  the  half-ellipsoid  and  Griffith 
notch.  An  estimate  of  the  amount  of  this  modification  can  be  made  for  the 
Griffith  notch  as  follows:  We  note  that  asa  + ® the  ellipsoid  approaches 

an  elliptic  cylinder,  so  that  (23)  should  be  valid  for  the  cylinder.  In  the 
case  of  the  half  cylinder  (Griffith  notch)  of  Figure  4b,  it  has  been  shown 
(Paris  and  Slh,  1965)  that  the  free  surface  causes  the  stress  concentration 
factor  for  tension  to  be  Increased  by  a factor  of  1.12,  while  that  for  shear 


max 


/ 

T(1  + V)  + Y 

*2 « - *>2  * (^y  J 

(24) 


This  equation  differs  appreciably  in  appearance  from  that  applying  to 
the  through-crack  (17a).  Among  other  things,  it  depends  on  Poisson's  ratio. 
However,  for  the  range  of  values  appropriate  to  glass,  v - 0.2  to  0.3,  the 
relative  contribution  of  the  applied  tensile  and  shear  stresses  to  the  maxi- 
mum tensile  stress  on  the  surface  of  the  cavity  is  almost  the  same. 

The  present  method  of  determining  failure  probability  curves  for  various 
stress  states  is  to  solve  the  integral  in  (6)  numerically  over  the  entire 
range  of  applied  stress.  The  data  for  equiblaxlal  tension  are  used  to 
determine  dN/docr  via  (12).  While  it  is  not  necessary  to  assume  a Welbull 
form  for  Pf(o),  several  Welbull  fits  to  the  date  were  investigated.  The  best 
fit,  however,  was  a piecewise  function  defined  over  three  regions  of  the  applied 
stress  range  (Figure  2)  and  used  exclusively  in  the  computations  of  the 
present  theory. 
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The  fraction  u>/ir  is  a function  of  the  applied  stress  state  and  o cr  and 
is  determined  in  the  following  manner.  The  maximum  tensile  stress  on  the  sur- 
face of  the  crack  (equations  17,  21  or  24)  can  be  written  as 


o « -r  o g(9,K) 
max  £ x ® 


(25) 


where  g(0,K)  is  the  function  of  crack  orientation  0 and  stress  state  K appro- 
priate to  the  crack  model  and  0^  is  the  applied  stress.  The  failure  criterion 

is  £o  2o  or  a g(0,K)  > o . The  fraction  w/n  is  the  portion  of  the 

max  cr  x cr 

range  0 < 0 < tt/2  for  which  g(0,K)  £ a Jo  , and  in  general  can  be  determined 

cr  x % 

from  the  roots  of  the  nonlinear  equation  g(0,K)  ■ o /o  at  each  value  of 

0 10  . 
cr'  x 

Figure  5 compares  the  failure  probability  curves  for  stress  ratio  1:0 
deduced  from  the  experimental  data  for  stress  ratio  1:1  using  (13),  (21)  and 
(24).  The  results  using  (21)  are  somewhat  closer  to  the  experimental  data, 
but  not  much.  A correction  of  unknown  magnitude  is  needed  to  account  for  the 
presence  of  the  free  surface.  Moreover,  it  is  doubtful  whether  a fracture 
criterion  should  be  based  on  a maximum  stress  occurring  at  the  end  of  the 
major  axis.  If  the  material  strength  is  exceeded  there,  the  crack  should 
lengthen  along  the  surface  but  not  penetrate  in  from  the  surface  in  the  manner 
required  to  partition  the  specimen.  Accordingly,  we  conclude  that  (24)  is 
the  more  appropriate  fracture  criterion. 

Use  of  the  critical  stress  concentration  factor  for  the  fracture  crite- 
rion leads  to  the  same  conclusion.  The  critical  stress  concentration  factors 


for  a notch  subjected  to  tension  and  out-of-plane  shear  are  (Paris  and  Sih, 
1965) . 


and 


K^j.  - s/irb  (27) 

respectively,  where  b is  the  appropriate  ellipsoid  axis  for  a Griffith  notch 
in  our  notation.  As  a result,  along  the  plane  of  symmetry  of  the  notch 

c, ' 1-n\Jh T <28»> 

T ‘y/ft  S ««>) 

(plane  strain)  (28c) 


a - vo 
x z 


Combining  these  to  get  the  maximum  principal  stress  as  a function  of  r by 
using  Mohr's  circle,  we  get 


max 


■ ia2v^  h1 + ^ + (ife)' 


(29) 


The  expressions  in  square  brackets  in  (24)  and  (29)  are  identical,  and  this 

is  all  that  is  used  to  find  the  relative  contribution  of  shear  to  crack  failure. 

It  is  of  interest,  however,  to  note  that  (24)  and  (29)  are  identical  when 
1 2 

r - - p - C /8b,  where  p is  the  minimum  radius  of  curvature  of  the  Griffith 
notch.  According  to  Mirandy  and  Paul  (1975),  this  is  the  value  of  r at  which 
the  solution  in  the  neighborhood  of  a closed  crack  tip  is  equal  to  the  value 
of  the  maximum  surface  stress  of  an  open  notch. 

We  find,  then,  that  agreement  with  experiment  is  improved  when  the  con- 
tribution of  shear  to  crack  fractures  is  taken  into  account,  but  that  this  is 
not  sufficient  to  bring  theory  into  really  good  agreement  with  experiment. 

We  oust  therefore  seek  other  reasons  for  the  remaining  discrepancy. 

One  possible  source  of  the  discrepancy  is  the  fact  that  the  data  for 
failures  at  a stress  ratio  1:1  cover  a somewhat  lower  stress  range  than 
failures  at  stress  ratio  1:0.  Thus,  an  analytical  expression  for  AN(acr) 
obtained  from  1:1  data  should  not  be  extrapolated  to  higher  values  of  o 


cr 
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However,  the  1:1  failure  stress  range  covers  most  of  the  1:0  failure  stress 
range  so  the  agreement  between  theory  and  experiment  in  this  large  overlap 
region  should  be  good,  but  it  is  not. 

Another  possible  source  for  the  discrepancy  is  the  assumption  in  the 
theory  that  crack  planes  are  always  normal  to  the  free  surface  of  the  material. 
This  explanation  also  fails.  It  has  been  shown  (by  methods  to  be  reported 
elsewhere)  that  the  1:0  failure  curve  deduced  from  1:1  data  is  the  same  using 
the  surface  crack  theory  described  here  as  it  is  using  the  volume  distributed 
crack  theory  of  Batdorf  and  Crose,  employing  a normal  stress  failure  criterion 
in  both  cases. 

The  final  possible  explanation  to  be  discussed  here  is  that  the  assumption 
of  uniform  distribution  of  crack  orientation  is  not  valid  for  the  pyrex  tubes 
tested.  There  is  no  a priori  method  of  predicting  what  form  an  anisotropic 
distribution  should  take.  A crude  but  useful  check  on  the  hypothesis  of 
anisotropy  is  to  investigate  the  consequences  of  assuming  that  cracks  are 
uniformly  distributed  through  a given  range  of  angles  and  are  absent  outside 
this  range.  For  instance,  one  might  assume  that  all  crack  planes  are  within 
radians  of  the  axis  of  the  pyrex  tubes,  thus  occupying  a fraction  R - 2^/tt 
of  the  total  available  angular  orientation.  Figures  6 and  7 show  the  results  of 
assuming  shear-sensitive  Griffith  notch  cracks  for  R - 0.6  and  R - 0.75.  It 
is  evident  that  the  anisotropy  assumption  greatly  improves  agreement  between 
theory  and  experiment.  Unfortunately,  the  tubes  were  tested  in  simple  ten- 
sion in  only  one  (the  circumferential)  direction,  so  that  a direct  test  of 
the  anisotropy  hypothesis  is  lacking. 


Figure  6.  Failure  Probability  Curvet  for  Uniaxial  Tantion  Calculated  Using  tha  B(1 :1)  Fit  to  tha 
Equibiaxiai  Data  for  QriffHb  Notch  Cracks  Assuming  a Simp  la  Anisotropy  of  Crack 
Or tantstion  in  Which  tha  Cracks  arc  Uniformly  Dbtrtoutad  Over  Only  a Certain  Fraction 
R of  tha  Possible  Angular  Orientations.  Curves  for  R - .6  and  R - .75  are  shown  in 
comparison  to  tha  data  points. 
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V.  CONCLUDING  DISCUSSION 


Most  statistical  fracture  theories  are  weakest  link  theories  that  attri- 
bute the  spread  In  fracture  stress  exhibited  by  nominally  Identical  specimens 
to  the  presence  of  invisible  flaws,  generally  believed  to  be  cracks.  Nearly 
all  such  theories  explicitly  or  Implicitly  assume  that  only  the  component  of 
stress  normal  to  a crack  plane  contributes  to  the  failure  of  the  crack.  This 
is  a convenient  approximation.  Thep calculations  In  this  paper  confirm  the 
expectation  that  the  errors  Involved  are  rather  small.  For  Instance,  In  the 
case  of  the  pyrex  glass  for  which  the  theory  was  evaluated  numerically,  the 
ratio  of  biaxial  to  uniaxial  stress  for  - 0.5  was  calculated  to  be  1.44 
for  shear-insensitive  cracks  and  1.38  for  shear-sensitive  cracks. 

Nevertheless,  since  the  assumption  of  shear  insensitivity  always  over- 
estimates the  biaxial  stress  failure  probability  calculated  from  uniaxial 
data,  an  improved  result  can  be  obtained  for  biaxial  stresses  in  the  tension- 
tension  quadrant  by  including  the  effects  of  shear.  The  present  paper  shows 
how  this  can  be  accomplished  in  the  case  of  surface  cracks.  A future  paper 
will  show  how  to  treat  the  case  of  volume-distributed  cracks. 
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